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Abstract 

We present numerical and analytical results for the thermodynamical properties 
of the spin-1/2 Heisenberg chain at arbitrary external magnetic field. Special em- 
phasis is placed on logarithmic corrections in the susceptibility and specific heat at 
very low temperatures (T/J = 10~ 24 ) and small fields. A longstanding controversy 
about the specific heat is resolved. At zero temperature the spin-Peierls exponent is 
calculated in dependence on the external magnetic field. This describes the energy 
response of the system to commensurate and incommensurate modulations of the 
lattice. The exponent for the spin gap in the incommensurate phase is given. 
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1 Introduction 



The seminal model for correlated quantum many-body systems is the well-known Heisen- 
berg model. In one spatial dimension the system with nearest neighbour exchange is 
exactly solvable for the spin-1/2 case [1]. The elementary excitations have been classified 
as spinons a long time ago [2,3]. Despite the integrability a comprehensive theoretical 
treatment is still missing. A notorious problem is posed by the calculation of correlation 
functions even at zero temperature, not to mention the corresponding properties at finite 
temperature. However, the asymptotics of the correlations in particular in the ground- 
state and at low temperatures are quite well understood by a combination of Bethe ansatz 
[4] and conformal field theory [5,6]. At exactly zero temperature the spin-1/2 Heisenberg 
chain shows algebraic decay in its correlation functions and thereby constitutes a quantum 
critical system. 

In recent years several experimental systems have been synthesized which realize very 
closely quasi one-dimensional isotropic Heisenberg chains [7,8]. In experimental studies 
the two-spinon continuum characteristic for the one-dimensional geometry has been veri- 
fied in the dynamical structure factor, see for instance [9] . Some other properties probing 
the critical properties of the system comprise corrections to the low-temperature asymp- 
totics in for instance the susceptibility [10]. Typically these corrections are logarithmic 
singularities [11,12,13,14] even without impurities. The origin of this behaviour can be 
traced back to the existence of marginal operators (in terms of statistical mechanics). 

Another aspect of the quantum criticality of the spin-1/2 chain is its instability towards 
a structural transition driven by long-ranged quantum fluctuations. Such instabilities 
generally known as Peierls or spin-Peierls transitions were originally treated in [15] for the 
case of free fermions (showing only marginal instability). In these days the understanding 
of the truly interacting case has increased in qualitative as well as quantitative aspects 
[16,17]. 

In Section 2 of this paper an approach to the thermodynamics on the basis of the 
physical excitations (spinons) is applied to the study of several thermodynamical prop- 
erties [18] (see also the related approach by [19,20,21,22,23,24,25,26]). The main results 
concern the numerical and analytical calculation of the universal leading and next-leading 
corrections in the susceptibility. In particular the resolution of a controversy between ex- 
act [27], field theoretical [28] and numerical approaches [12,13,14,29] about the specific 
heat at low temperatures is presented which is absolutely free of fitting procedures in 
the analysis of the zero-temperature limit. Our approach permits the numerical study of 
extremely low-temperatures down to T/J = 10~ 24 . Furthermore the analytic treatment 
reveals the origin of the logarithmic corrections in certain properties of the spinon-spinon 
scattering phase. In Section 3 we calculate the critical exponents of correlation functions 
at zero temperature in the presence of an external magnetic field. On the basis of scaling 
relations the spin-Peierls exponents describing the gain in magnetic energy and the open- 
ing of gaps (even in the incommensurate phase) are derived. This is expected to be the 
basis of a more accurate treatment of the temperature-field phase diagram of spin-Peierls 
systems. 
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2 Thermodynamics and logarithmic corrections 



We investigate the thermodynamical properties of the antiferromagnetic isotropic spin-1/2 
Heisenberg chain 

L 

H=2JJ2$A+v (1) 

j=i 

We note that the elementary excitations ( "spinons" ) have quasi-linear energy-momentum 
dispersion with velocity 

v = n,J. (2) 

The free energy per lattice site of the system is given by the following set of non-linear 
integral equations [18] for auxiliary functions o, and 21 = 1 + a 

logoff) = V <t> + / \k(x — y) log2l(?/) — k(x — y — in + ie) log2l(y)l dy, (3) 

coshx 

and the corresponding equation for a, and 21 = 1 + a is obtained from (|3|) by exchanging 
a <-> o, 21 <->• 21, and i — > —i, <fi —> — 0. Finally, <j) is given by the external magnetic field h 
through <p = Bh/2. The integration kernel is defined by the Fourier integral 

k{x)= H ^k cos{kx)dk - (4) 

In terms of the solution 21 and 21 to the integral equations the free energy is given by 

i rM^. (B) 

Z7T J _ QO cosnx 

These equations are readily solved numerically for arbitrary fields and temperatures by 
utilizing the difference type of the integral kernel and Fast Fourier Transform. In Fig. 
the specific heat c(T), susceptibility x(T) and correlation length £(T) (computed on the 
basis of [18]) at zero field are presented. Note the divergence of £(T) at low temperatures. 
In Figs. and |3| the results for the specific heat and susceptibility are shown for various 
magnetic fields below as well as above the saturation field h c = 4J. Some qualitative 
aspects of the curves in dependence on the fixed external field can be understood in the 
picture of spinon excitations. The external magnetic field acts very much like a chemical 
potential for the spinons for which there are particle and hole like excitations. At zero 
field the bands of the particle and the hole type excitations are identical, however for 
sufficiently strong field h the band widths are considerably different resulting into two 
maxima in the specific heat at different temperatures. 

In Fig. | (a) the susceptibility x(T) is plotted down to much lower temperatures show- 
ing the onset of an infinite slope at T = as noticed in [11]. In Fig. ^ (b) a low tempera- 
ture analysis of xCO down to T/J = 10~ 24 is presented giving numerical evidence about 
leading and subleading logarithmic corrections in temperature. Here we like to present a 
completely analytical low-temperature analysis (for h «T) and defer comments about 
the numerical analysis to later. 
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Figure 1: Plots of the specific heat c(T), susceptibility xCO ( m units of 1/J) and corre- 
lation length £(T) at zero external field versus T (in units of J). 



Analytical study of low-temperature asymptotics 

Because of the driving term on the RHS of (§) it is convenient to introduce the following 
scaling functions 



a £ (x) = a(x + £), £ = log(vr/3), 
0£(x) = o(— x — £), 



(6) 



approaching well-defined non-trivial limiting functions in the low-temperature limit which 
satisfy 

loga £ (x) = -2e- x + O(l//3) + + <Mx) 

/oo 
[k(x - y) log Slfi(y) - k(x - y - in + it) log 2t £ (y)] dy, (7) 

where i/j£,(x) accounts for the contribution of the functions 21 and 21 on the negative real 
axis 



k(x + y + 2£) log 2l £ (y) - fc(ar + y + 2£ - ivr + ze) log 2l £ (y) 



dy. (8) 



Using these functions the free energy at low temperatures takes the form 

1 



/ = e - 



e x log 



2l £ 2l £ 2l £ 2l £ (x) dx + 0(l//3 3 ) 



(9) 
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Figure 2: Depiction of the specific heat of the Heisenberg chain for different values of 
the magnetic field (a) h/J = 0, 0.5, 2.0, (b) h/J = 2.5,. ..,3.5, (c) h/J = 4.0,. ..,9.5. 
Note the linear T dependence at low temperatures T for fields less than the saturation 
value h c j J = 4, see (a) and (b). For h = h c the specific heat is proportional to T 1 / 2 for 
sufficiently low T, and for even stronger fields a thermodynamically activated behaviour 
develops, see (c). 



In order to analyse the low-temperature asymptotics in detail we perform the following 
manipulation 

/oo 
[(log ^(log a*) - (logo*) (log St*)'] dx + ! conj\ 
£ 

/oo 
[(-2e~ x + + ^(x))'(log2l £ ) - (-2e" x + + ^(x))(log2l £ )'] dx + 'conj', 
-£ 

/oo 
(2e" x + ^(x)')(log2l £ )dx - 01og (1 + e 2< ^) + 'conj', (10) 
-£ 

where in the second line we have inserted (0) and in the third line we have performed 
an integration by parts with an explicit evaluation of the contribution by the integral 
terminals using the limits 

loga £ (cx)) = +20, log2l £ (oo) = log(l + e +2 ^), 

loga £ (oo) = -20, log2t £ (oo) = log(l + e- 2 ^). (11) 
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Figure 3: Depiction of the susceptibility (in units of 1/J) of the Heisenberg chain for 
different values of the magnetic field (a) h/J = 0, 0.5, 2.0, (b) h/J = 2.5, ...,3.5, (c) 
h/J — 4.0, ...,9.5. Note the finite value at zero temperature for fields less than the satu- 
ration value hj J = 4, see (a) and (b). For h = h c the susceptibility diverges at T = 0, 
and for fields larger than h c the susceptibility drops to zero, see (c). 



Next, we show that the LHS of © can be evaluated explicitly as it is a definite integral 
of dilogarithmic type with known terminals 



LHS = 2/ (logo £ )'log2l £ dx-loga £ (log2l £ )' 
J-s. 

= 2 / log(l + e z )dz-201og(l + e 20 ) + 2 

J — oo 
7T 2 

3 ' 



+ 'conj 



-2d> 



log(l + e z )dz + 20 log(l + e" 2 ^) 



'121 



and the final result for the definite integral in the second line was obtained by standard 
techniques (it is independent of (ft and the particular choice = "directly" leads to the 
final result). 

Combining (|T0"D and ([Tj^) we arrive at a formula particularly suited for studying the 
asymptotic behaviour of (|S|) 



7T 



I / e- x log[2l £ 2l £ ]dx = y + 



2 / [^(x)' log 2l £ + VbW log Sa c ] dx. (13) 



From the general expression of the free energy / (|5|) in terms of the functions 21 and 
21 the entire functional dependence of 21 and 21 matters even for the calculation of the 
asymptotic behaviour of / at low temperature. The usefulness of (|1^) lies in the fact that 
here only the asymptotic behaviour of 21 and 21 at large spectral parameter enters. The 
reason is simple to understand: for any finite argument x the contribution is of order l/£ 2 
due to the i/)q(x) factors. Here we see that the low-temperature asymptotics (and further 
below the occurrence of logarithmic corrections) are intimately linked to the algebraic 
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Figure 4: Plots of (a) zero-field susceptibility xCO versus T/J showing a logarithmic 
singularity at T = and (b) the next-leading logarithmic correction (x/xo — 1 — l/2£)£ 2 
versus log£ (~ log log J/T) down to T/J = 10~ 24 . Note the linear dependence on log£ 
for low temperatures. 



decay of the integration kernel k(x) 

*<*> = OT' (14) 

We like to add that the kernel is related to the S'-matrix of the elementary spinons [3] by 
k(x) = £\ogS(x). 

Asymptotics of the susceptibility 

In a first step we determine the asymptotics of the auxiliary functions from (|7]) where 
we drop integrations over negative arguments and for large arguments we use ci£ ~ l/ci£ 
(hence logSl^ — logSle = logae) and furthermore we equate the kernels &(...) with argu- 
ments differing by iri leading to 



log az(x) = <p + ipz(x) + k(x- y) log a £ (y)dy, 

Jo 

POO CI POO 

= + k{x-y)log— — -dy+ / k(x - y) logo £ (oo)dy, 
Jo a £ (oo) J 

= <p + ipz(x) + log a& , X * [ k(x - y)dy 

Jo 



k(x - y)loga £ (oo)dy - / k(x-y) log az(oo)dy, 



oo J — oo 



+ ^&( x ) + iloga £ (x) - J 



log a&(y) i 
A{x-yf ay ' 



+ ?Pz(x) + ^loga £ (x) - — . (15) 



In the first integral of the second line we have used the fact that dominant contributions 
occur for the integration variable y close to x. In the subsequent calculations the asymp- 
totic behaviour of the kernel and also the total integral J k(x)dx = 1/2 were used. Finally, 
employing the limits ([TT|) we arrive at 

loga £ (a;) = 2<f> - t (16) 
x 

In the second step we insert (H) into the integral on the RHS of (0), then we insert 
([[ID and calculate the resulting integral by dropping terms leading to contributions of 
order 0(l/£ 2 ) and higher 



-£ 



[ip £ (x)' log Ql 2 + ^ z {x)' log 2l £ ] dx 



= +#2 r r k>{x +y+2 ^( i - Yx) i 1 - Y y ) ^ 

= [ k(x + 22) (l - I) dx = -0 2 /°° ( X _ I) dx 

Hence we obtain the dependence of the free energy on the magnetic field at low temper- 
ature 

4 £ e- >og[ Wx = £ + 2*» (l + -L _ ) , (18) 

implying logarithmic corrections in T for the free energy and the zero-field susceptibility 

/ = e -^-^(l + ±- l -^§-), £ = log(7r/T), (h«n (19) 
17 6w Ativ V 2£ (2£) 2 y ' b\ / j, \ 

1 ( 1 log(2£)\ . . 

X = w-. U+^-^r ■ (20) 



2tto V 2£ (2£) 2 y 

This result is in complete agreement with the numerical analysis in Fig. |] (b) and with 
[30,29] (if the system size L is replaced by the inverse temperature 1/T). The (non- 
universal) temperature scale To entering the logarithms of the correction terms will be 
studied elsewhere [31]. 



Asymptotics of the specific heat 

Next we calculate the higher order low-temperature asymptotics of the Heisenberg chain 
in zero external magnetic field. Here the functions ci£ and cic are related by complex 
conjugation. We therefore write 

log2t £ (x) = logf £ (x) = R £ (x) +il a (x), 

log2t £ (x) = loga £ (z) = ifc(x)-iJ £ (z). (21) 
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Again we insert (|8]) into the integral on the RHS of ([13]), then we insert ( |21| ) and expand 
k(x + 7rz) = k(x) + k'(x)iri — k"(x)7T 2 /2 obtaining 



[Mx)' log SM-jc + VsO)' log 2l £ ] cfe 



oo ^oo 



-£ ./-£ 



OO /"OO 



dxdy[7r 2 k"'(x + y + 2Z)Rz(x)Rz(y) + 4A/(x + y + 2£)/ £ (x)/ £ (y) 

- 2ttA;"(x + y + 2£) (R & (x)I & (y) + I a (x)R £ (y)) ] , 
cfody [vr 2 A;'(x + y + 2£,)R' s ,(x)R' z {y) + Ak'{x + y + 2£)I £ (x)I £ (y) 

+ 2nk'(x + y + 2£) (R' 2 (x)I £ (y) + h{x)R'M) ] , (22) 



where in the last step we have performed an integration by parts. Next, the integral is 
dominated by contributions with finite integration variable x and y for which k'(...) can 
be replaced by a constant yielding 



-£ J-s, 



RHS 



OO POO 



16£ 3 



dxdyfc&six) + 2I z {x)] [nR'M + 2/ £ (y)] 
vr.R + 2/) 2 , 



16£ 3 j_ £ 
1 



(23) 



where we have used the abbreviations 

R = 



dxR'^x) = log 2, 



/OO 
dxl^(x) — > numerical evaluation, 
-£ 



(24) 



where i? can be calculated easily from the asymptotics of R^{x), I however has to be 
evaluated by numerical calculations. For the free energy we obtain 



I / e- x log[2l £ 2l £ ]dx = ^ 

-£ 6 



(25) 



where (i? + -Jj is found to be 1 ± 10 6 . This leads to 

3 1 



7T 



(26) 



confirming the field theoretical prediction in [28] about the amplitude of the l/(logT) 3 
correction which was argued to be universal however in disagreement with Bethe ansatz 
calculations in [27]. The reason for the "failure" of the treatment in [27] (instead of 
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Figure 5: Plots of the estimators of the leading logarithmic corrections (a) for the suscep- 
tibility (x/xo ~ 1)£ an d (b) the specific heat [3vc(T) / 'irT — 1]|£ 3 versus log 10 (T/J). Note 
the non- monotonous behaviour of the higher order terms in (b). 



3/8 the value 0.3433 was obtained) is subtle and lies in an inappropriate use of the 
Euler-MacLaurin formula which unfortunately is common practice in the treatment of 
critical lattice systems. The numerical treatments of the Bethe ansatz equations up to 
iV = 16384 [12,13,14] were still plagued by higher order logarithmic contributions leading 
to erroneous conclusions with respect to the analysis of the specific heat. For an illustration 
of these higher order terms see Fig. |5|. For the case of the specific heat an appropriate 
fit of finite size data seems very difficult. We note, that the careful numerical treatment 
[29] of finite systems up N = 16384 employed an appropriate fit algorithm taking into 
account the proper higher correction terms. Our results agree with the conclusions of 
[29]. We emphasize that our treatment is completely free of any fit procedures. The type 
of logarithmic singularities has been derived analytically and the amplitudes have been 
calculated (in the worst case) numerically which did not face any extrapolation problems. 



3 Groundstate properties of homogeneous and 
dimerized S=l/2 Heisenberg chains 

Next we treat the critical groundstate properties of the homogeneous Heisenberg chain 
by means of an application of the exact Bethe ansatz solution and conformal field theory, 
see [32] and references therein. 
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Bulk properties 

Directly from the Bethe ansatz equations for spin rapidities the (linear) integral equation 
for the groundstate density function p{v) is derived [4,32] 

p{x) = p (x) - k(x- y)p(y)dy, (27) 



J -K 

where the so-called bare density po and the kernel k(x) are defined by 

1 2 , , , 1 4 

Po = 5 , k(x) = = -. 28 

FU 2vr x 2 + 1 v ; 2vr x 2 + 4 v ; 

Here we have as yet left unspecified the magnetization m of the state (or equivalently the 
external magnetic field h) corresponding to some integration terminal K with subsidiary 
condition 

1 f K 

m = - - / p{x)dx, 2k F = (1 - 2m)vr, (29) 



2 



where (twice) the Fermi momentum kp is the momentum at which elementary spin exci- 
tations become soft. In terms of the density function the groundstate energy e takes the 
form 

-K 

€o{x)p(x)dx, (30) 

-K 

where e (x) is the bare energy of spin excitations 

4J 

e °( x ) = TT^^ = X » ( 31 ) 

x A + 1 

and the bare charge function ^o( s ) has been introduced for later purposes. 

In the same way as (|27j ) defines a function p for the inhomogeneity poj h gives rise 
to the dressed energy e(x) and the dressed charge £(x) related to (|31~1). (Note that the 
auxiliary function a(x) of the thermodynamic treatment is related to e(x) in the low- 
temperature limit via log a = (3e and (§) turns into the linear integral equation for the 
dressed energy.) As the integration operator in ( p7|) is selfadjoint we have an alternative 
relation to ( p0|) and (|29|) for the groundstate energy e and magnetization m 

k ^ f K 



e(x)p (x)dx, m = / £(x)p (x)dx. (32) 

k 2 J_ K 

For relating the terminal K to the magnetic field h we note that the bare energy 
in the presence of an external field is given by e (x) + h = e (x) + h^ (x). Hence the 
corresponding dressed energy is e(x, h) = e(x) + h^(x). A general property of the dressed 
energy is that it has to vanish at K yielding 

e{K,h) = 0*h = -^l. (33) 
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For any value of K the functions e(x), £(x) can be evaluated (at least numerically) from 
an integral equation of type (^7|), then the magnetic field h, and the corresponding mag- 
netization m and groundstate energy e — mh can be calculated utilizing fl33|) and ([32|) . 

Correlation functions 

We are particularly interested in the determination of scaling dimensions (critical expo- 
nents) for arbitrary external fields. For instance the spin-spin and energy-energy correla- 
tions decay algebraically with exponents x s and x e 

(- 1 ) r / v cos(2A; F r) C 
(S S r ) ~ e • e r ) ~ + -, (34 

where the energy operator is e r = S r S r +i. In general the exponents are directly given in 
terms of the dressed charge function at the terminal K (see [32]) 

x = ^(KfS 2 + ^Kfd^ (35) 

where S and d are integers corresponding to the quantum numbers of the particular 
operator (occurring in the correlation function); S equals the magnetization and d equals 
the lattice momentum in units of 2k p. For spin and energy operators (correlations) we 
have (S,d) = (1,0) and (S,d) = (0, 1), respectively Therefore we obtain in terms of the 
dressed charge 

x ' = 4£(Ky>> X * = W. (36) 



We want to note that ( |35| ) is obtained from predictions by conformal field theory 
about finite-size effects for the spectrum of the system on large lattices L with periodic 
boundary conditions 

E x -E = 2ttj(x + 8(L)), (37) 

where E Q is the groundstate energy, E x the low-lying energy level corresponding to the 
operator with dimension x, v the velocity of the elementary excitations, and S(L) is a 
function with zero limit for L — > oo. The evaluation of (|37|) can be carried out analyt- 
ically yielding (|35|) (see for instance [32] and references therein). Of course, the dressed 
charge £ has to be calculated numerically. Alternatively, similar results can be obtained 
by numerical calculations on the basis of the Bethe ansatz equations for finite systems and 
(B7p, however involving an explicit extrapolation analysis with the usual disadvantages, 
see [33] and for finite external fields [34]. 
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Perturbations of the pure system 



The algebraic decay of the correlation functions in the groundstate of the spin-1/2 Heisen- 
berg chain constitutes quantum criticality. Starting from this observation we are going 
to argue that the response of the system to perturbations can be understood and evalu- 
ated in a renormalization group framework and in general shows non-integer exponents 
(in contrast to simple perturbation treatments which of course are not legitime). We 
want to examine the critical exponents of the Spin-Peierls transition which is a structural 
transition of the underlying lattice with a modification of the local exchange couplings 
proportional to the local lattice distortion 



H(S) = H + 



SjS j+ i 



H + 5H. 



(3* 



An application of scaling relations will provide a simple tool to understand some of the 
essential aspects of the theory of [16,17]. We remind the reader of the response of the 
free energy / and the correlation length £ of a critical, classical system in d dimensions 
perturbed by a relevant operator with RG eigenvalue y(> 0) 



classical system 
d dimension 



H{5) = H + 5-H' 
free energy Af <■ 
corr. length 



d 
5v 

i 



A quantum critical system in d dimensions behaves very much the same as it is equiv- 
alent to a d + z dimensional classical system, where z is the dynamical critical exponent. 



quantum critical system 

d dimension 

g.s. energy e, gap mo 



classical system 
d + z dimension 

/, r 1 



For conformally invariant chains we have d = z = 1. Furthermore the RG eigenvalue 
y is related to the scaling dimension x of the particular operator by 



x + y = d + z = 2. 



(39) 



Now equating the above relations we obtain for the response of the quantum chain in the 
ground state energy per site e and the excitation gap m (= at the unperturbed point): 



Ae 



2 

-S2-X. 



m 



l 

(J 2-*. 



(40) 



For the spin-Peierls transition we have to use x = x e = ^(K) 2 . In Fig. |6| the value of the 
groundstate exponent a = 2/(2 — x e ) is depicted and shows a variation from the value 
4/3 for zero external field to 2 at the saturation field. As the exponent is less than 2 
the gain in magnetic energy always wins in competition with the loss in elastic energy 
which is of order 5 2 resulting in a structural instability, cp. Fig. |7j Of course this result 
is completely consistent with [16,17] where bosonization results were used for the energy- 
energy correlations of the spin-1/2 XX Z chain in vanishing external field in combination 
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Figure 6: Depiction of scaling dimensions x s , x e , groundstate energy exponent a, and 
Fermi momentum 



with scaling relations for the polarizability. In our treatment we have correctly accounted 
for the field dependence of the exponents which is new. 

Our study is partly motivated by the general interest in a deeper understanding of the 
physics of Spin-Peierls systems. Here we have shown how a combination of established 
results by conformal field theory, RG theory and integrable systems yields the exact critical 
exponents. At the same time we intend to use our results for a more precise calculation of 
for instance the phase diagram of Spin-Peierls systems in dependence on temperature and 
magnetic field. Note that there are quantitative discrepancies in for instance the phase 
boundary between the incommensurate and the uniform phase as obtained experimentally 
[35,36] and theoretically [17], cf. the remark in [36]. In addition, we note that the occurence 
of a gap in the magnetic system (for fixed modulation of the exchange couplings) is a very 
natural consequence of our approach. Experimentally it is still an open question if a spin 
gap in the incommensurate phase exists. The T 3 contributions to the specific heat at low 
temperatures [37,38] are explained in terms of gap less "phasons" [39]. 

Up to now we have not specified any particular dependence of the perturbation field Sj 
on the site index j. The above argument (based on a relevant perturbation) is applicable 
just for the case of a simple algebraic decay of ({8i€i){5jej)) for (i — j) — > oo without 
any oscillations which imposes a matching condition: (i) the modulation is sinusoidal 
8j = 5cos(qj), and (ii) the wave number q is identical to that of the energy-energy 
correlation q = 2kp. (For sinusoidal modulations with different wave number the response 
in the groundstate energy has the exponent 2.) 

As is indicated by the results in Fig. |6| the external magnetic field leads to a filling of the 
spinon-Fermi sea and a change of exponents x = x(h) and Fermi momentum kp = kp(h). 
Consequentially an increase in h amounts to (i) a suppression of the SP-instability and 
(ii) an incommensurate lattice distortion (at least for sufficiently strong fields depending 
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total 



5 



Figure 7: Competition of gain in magnetic energy with loss in elastic energy E tota i(5) = 

+ ^elast = ~Cl5 a + C 2 5 2 . 



on the balance of energy gain in the commensurate versus incommensurate case). These 
qualitative features are observed in [35,36]. Let us consider in some detail the response 
of the magnetic system (for small magnetization m) to a sinusoidal modulation of the 
exchange couplings (with wavevector 2kp). We are interested in the effects of Umklapp 
processes which favour commensurate modulations. For m = (2kp = 7r) we expect 
£ma g = ~ciS a and for m ^ (2k F ^ vr) only £ mag = -2 Cl (5/2) a = -2 l - a Cl 5 a . Hence, 
the groundstate energy of the total system as a function of the magnetization m shows a 
jump at m = 0. This is the origin of a 1st order transition from the commensurate to the 
incommensurate phase at a finite critical field h c . Whether higher harmonics ("solitonic 
modulations") can actually change the order of the transition from 1 to 2 cannot be 
answered decisively in this approach. 

Lastly, we want to comment on logarithmic corrections to the above results. These 
corrections typically appear if marginal operators exist, i.e. if there are interaction terms 
compatible with the usual symmetries of the lattice system (lattice translation, spin rota- 
tion), not however with strict conformal invariance in the continuum limit. Furthermore, 
these interactions have to possess the scaling dimension x m = 2. From the required sym- 
metry properties only the case (S, d) = (0, 2) in fl35|) is singled out with scaling dimension 
A^(K) 2 = 4x e . With a glance to Fig.|6] we directly see that logarithmic corrections only 
occur for vanishing external field. For completeness we want to note the physical scenario 
at this point. The function 5(L) in (|37D , which typically decays algebraically now has a 
leading logarithmic contribution S(L) ~ 6/logL. The amplitude b also determines the 
multiplicative logarithmic corrections [28] to 
(i) the two-point correlation function of an operator o(r) 



(ii) the response of the critical system to a (relevant) perturbation with H' = Y^ r =i °( r ) 
opening a gap 



where in both cases x is the scaling dimension of the operator o(r) and 5 is the (small) 




(41) 




(42) 
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coupling parameter. In our case and vanishing external field we have x — 1/2 and 6 = 3/2 
[28] such that m ~ 5 2 / 3 /(log 5) 1 ' 2 . 

4 Conclusion 

In this work two aspects of the critical spin-1/2 Heisenberg chain were treated. First, 
the logarithmic corrections to the leading low-temperature properties in the specific heat 
and susceptibility have been presented. For this purpose we have developed a system- 
atic expansion of the non-linear integral equation governing the thermodynamics of the 
system. In physical terms, the low-temperature singularities are caused by the algebraic 
asymptotics of the S'-matrix of elementary excitations in dependence on the difference 
in the spectral parameters. This is consistent with the absence of logarithmic corrections 
for models with long-range interactions [40]. In addition, numerical results have been 
presented for the Heisenberg chain down to very low temperatures. 

Second, at zero temperature we have treated the perturbation of the Heisenberg chain 
by inhomogeneous couplings. Quite generally this problem can be solved by use of scaling 
relations. In our case the Spin-Peierls exponents in presence of an external magnetic 
field have been determined. We like to note that a similar reasoning is applicable to 
many questions about critical quantum systems perturbed by additional couplings, e.g. 
interchain couplings. 

In future work we want to combine the two approaches, calculation of finite temper- 
ature correlation lengths and exact scaling dimensions, to a more detailed description 
of the spin-Peierls scenario of the Heisenberg system, e.g the critical lines of the phase 
diagram. 
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